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(57) Abstract 

Raw image data is acquired as it is being produced by a magnetic resonance unit during a functional magnetic resonance imaging 
(fMRI) examination. Reconstruction, registration, and statistical analysis of the acquired raw image data are performed, and the results arc 
displayed immediately after completion of the acquisition of the raw image data. A parallel processing technique is used to achieve the 
real-time interactive fMRI. 
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REAL-TIME INTERACTIVE FUNCTIONAL 
MAGNETIC RESONANCE IMAGING 

FIELD OF THE INVENTION 

5 The invention relates generally to the real-time display of images for interactive 

functional magnetic resonance imaging (fMRI). 
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25 

BACKGROUND OF THE INVENTION 

fMRJ facilitates the study of dynamic physiologic processes by measuring the 

changes in magnetic resonance (MR) signal intensity during rapidly acquired serial images. 

fMRI studies, for example, can be used to determine the following: 1) relative cerebral 

30 blood volume (rCBV) and relative cerebral blood flow (rCBF) using dynamic contrast 

enhanced MRI (1-1 1); 2) absolute cerebral blood flow (CBF) using spin tagging techniques 

(12-16); and 3) localized increases of neuronal activity in response to neurophysiological 
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stimulation using blood oxygenation level dependent (BOLD) contrast (17-25). By using 
rapid scanning techniques, such as single shot echo planar imaging (EPI) or spiral imaging, 
each fMRJ study usually acquires large raw data sets, which are comprised of hundreds or 
thousands of images, in relatively short periods of time (18, 24, 25). 
5 These fMRI data sets are conventionally not reconstructed on the MR unit because 

reconstruction would interfere with or delay the start of subsequent fMRI experiments. 
Instead, the raw data is conventionally stored and subsequently transferred to an off-line 
workstation for analysis, which usually occurs long after the subject has left the MR unit 
(21-25). 

10 Typically, analysis of fMRI data sets usually take several hours to days (21, 25) 

before they are reviewed because of the delays in off-line reconstruction, registration, and 
subsequent statistical analysis of the raw image data. Time is also required to archive all 
unprocessed raw image data files and images. Although the most important part of this 
computational process lies in the ultimate analysis of the fMRI images, for which various 

15 software packages are available (33, 34, 35, 36), the bulk of the initial time commitment is 
spent on the reconstruction, registration, and storage of the images. As a result of this 
conventional process, even a quick review of activation maps usually occurs long after the 
subject has left the MR scanner. Due to delays imposed by the MR unit, researchers have 
tended to limit the number of experimental trials acquired during a one-hour scanning 

20 session, in part because the researchers are unsure of the reliability and integrity of their 
studies. 

Moreover, the fMRI results are often inadequate due to motion artifacts or poor 
compliance because of inadequate understanding of the paradigm. To obtain adequate 
fMRI results when an initial fMRI exam has gone awry, the subject needs to reschedule 
25 another fMRI exam, return to the fMRI facility, undergo the additional fMRI exam, await 
the results of the additional fMRI exam, and possibly reschedule yet another fMRI exam if 
the results are again poor. Because of the significant costs of MRI exam time, the 
increased backlog of MRI exams, and the increased discomfort and displeasure of the 
subject, this conventional approach of repetitive studies is unsatisfactory for both the 
30 researcher and the subject. 

Real-time fMRI approaches have been proposed previously. For example, Cox et 

al. (21) proposed an algorithm that recursively computes the correlation coefficient of an 

.5. 
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image sequence with a reference vector and simultaneously removes any undesired linear 
trend from the data. Activation maps were displayed using a predetermined statistical 
threshold. This process was performed at various time points during a BOLD EPI study 
acquired at a single location while the subject performed a finger tapping experiment. The 

5 authors indicated that their algorithm was not a complete substitute for post-processing, 
where multiple statistical tests can be applied to the raw data sets for detection of the 
activation and suppression of the artifacts. However, the authors indicated that the benefit 
of their approach is to look primarily for spuriously correlated artifacts. 

With the approach of Cox et al. (21), the temptation exists to truncate an 

10 experiment early because the statistical maps generated appear to give the results that the 
investigators required. This temptation should be avoided. A decision to halt a study early 
could result in an arbitrary number of epochs collected, which would confound any type of 
intra-subject or inter-subject comparisons. Moreover, any interim analysis of an fMRI 
study as a consequence of a real-time display should impose a more stringent (higher) 

15 threshold level for significance to define activation. 

As another example of prior art attempts, real-time reconstruction and display of 
MR data have been used successfully to direct single slice MR fluoroscopic based imaging 
techniques to localize anatomical structures of guide and monitor intervention procedures 
(25, 39-41). Further, Kerr et al. (26) demonstrated the utility of a real time interactive MRI 

20 on a conventional scanner with abdominal and cardiac imaging. 

SUMMARY OF THE INVENTION 

With the invention, an fMRI examiner can immediately evaluate the results from an 
fMRI examination of a subject while the subject is still in the MR unit. If the results are 
25 poor, the examination can be performed again without the need or expense of having to 
reschedule the subject for another fMRI examination; Further, if the results dictate 
performing additional fMRI examinations, these examinations can be immediately 
performed without the need or expense of having to reschedule the subject for another 
fMRI examination. 

30 For example, with the invention, results of a whole brain BOLD activation 

experiment with isotropic voxels (17, 18, 24) could guide the planning of a subsequent 
perfusion experiment (i.e., dynamic contrast enhancement or spin tagging). These types of 
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studies, all performed within a single fMRI exam session, allow the examiner to review the 
BOLD activation map and then to either vary the work load of the task or repeat the study 
multiple times to observe changes associated with reproducibility, learning, or precision. 
The advantages of real-time interactive fMRI include: 1) examining the data for 

5 artifacts that may contaminate the images (i.e., physiological motion); 2) developing or 
modifying activation paradigms quickly, allowing the examiner to test an hypothesis while 
the subject is still in the MR unit; and 3) developing new paradigms on-line, and thereby 
making fMRI a more flexible neurological tool. Furthermore, the immediate display of the 
analyzed results of BOLD, perfusion (i.e., spin tagging) or dynamic contrast fMRI studies 

10 also allows for the evaluation of the subject to a provocative pharmacological challenge 
(e.g., the Wada or sodium ambitol test). Real time analysis of an fMRI exam also permits 
a direct comparison to the baseline condition of the subject, and may provide valuable 
information relating to neurological integrity and functional status of the subject. With 
interactive fMRI, sequential fMRI exams can be designed in which, for example, the 

15 results of an initial BOLD activation study guides other functional (i.e., spin tagging, 

dynamic contrast) or metabolic imaging (37) examinations to a specific area of interest or 
to repeat the paradigm with changes, such as higher spatial and/or temporal resolution (38). 

Interactive fMRI can also be useful to dynamic contrast studies. In particular, 
dynamic contrast enhanced fMRI studies are usually performed once during an imaging 

20 session and are being used as part of an early evaluation of an acute ischemic event or an 
evaluation of central nervous system (CNS) pathology (1,6). These results may be useful 
in monitoring the treatment effects of anti thrombolytic agents in an acute stroke (4), and 
have been shown to be useful in grading primary central nervous system (CNS) 
malignancies (16). With the invention, analysis of the calculated rCBV, rCBF, and time- 

25 to-peak (TTP) maps after completion of the acquisition of the raw image data may provide 
further insight into the efficacy of novel therapies. 

The success of an interactive fMRI based physiological interview depends on the 
investigator being able to: 1) review processed and statistically analyzed image maps 
shortly after the completion of the fMRI acquisition; 2) determine if the study needs to be 

30 repeated; or 3) if a modification is required to the paradigm, further elucidate cerebral 
response. This sequence of events should occur while the subject is in the MR unit in 
order to take advantage of the experimental conditions. With the invention, the 
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neuroscientist is provided with an analysis of large fMRI data sets within seconds after 
completion of the MR acquisition. 

The invention, in its preferred embodiment, costs less than most major MRI 
hardware upgrades, while simultaneously providing researchers with a more efficient use 
5 of MR scanner time, and ultimately increasing the number of clinical and research studies 
performed. 

Moreover, the invention facilitates the on-line planning and evaluation of various 
protocols (e.g., spectroscopic imaging, dynamic contrast CBV, arrival time maps, and 
perfusion CBF measurements. The invention can be used in conjunction with other 

10 modalities to study the effect of carotid stenosis on cerebral perfusion ( 1 1 ), or become part 
of a work-up for a brain attack (1,5, 42, 43). Further, the invention can be used for 
adaptive task tuning (e.g., increasing the stress-load) or repeating unsuccessful scans (e.g., 
subject motion) used for brain mapping (12, 18), or pre-surgical planning of eloquent 
cortex relationship to a surgical target (3, 25, 23, 44), and other therapeutic interventions, 

15 as will become apparent to one of ordinary skill in the art. Moreover, the invention can be 
used to integrate fMRI techniques into clinical practice, and as a valuable tool for 
improving the care of patients with neurological and psychiatric disorders. 

It is an object of the invention to provide a parallel processor hardware system, 
along with post-processing techniques, to display statistical results of an fMRI exam co- 

20 registered on anatomic images and to allow for the rapid evaluation and preliminary 
interpretation of the fMRI exam. 

It is another object of the invention to provide a hardware and software 
configuration using parallel processors that cost less than most MR hardware upgrades and 
that rapidly performs the following tasks: 1) performs MR acquisition by capturing the 

25 raw image data on-the-fly; 2) reconstructs the image files; 3) registers all volumes; 4) 
performs a predefined statistical analysis for either a BOLD fMRI study or a dynamic 
contrast study; and 5) displays relevant maps within seconds after the completion of the 
MR acquisition. 

The above objects and advantages are achieved by the invention, which includes a 
30 method, an apparatus, and an article of manufacture for the real-time display of the results 
from an interactive fMRI examination. 

-8- 
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In particular, the invention includes a method for outputting results from a magnetic 
resonance image examination. The method comprises: obtaining first and second 
acquisition volumes from a magnetic resonance image scanner during the magnetic 
resonance image examination; processing the acquisition volumes from the magnetic 

5 resonance image scanner to obtain the examination results; and outputting the examination 
results! The processing of the acquisition volumes comprises: reconstructing the first 
acquisition volume using a first processor to obtain a first reconstructed volume; 
reconstructing the second acquisition volume using a second processor to obtain a second 
reconstructed volume; and registering the second reconstructed volume using a third 

10 processor. 

Alternatively, in the method of the invention, the processing of the acquisition 
volumes comprises: reconstructing the first acquisition volume using a first processor to 
obtain a first reconstructed volume; reconstructing the second acquisition volume using the 
first processor to obtain a second reconstructed volume; and registering the second 

15 reconstructed volume using a second processor. 

Further, the invention includes an apparatus for processing data from a magnetic 
resonance image examination. The apparatus comprises: a magnetic resonance image 
scanner for producing a plurality of acquisition volumes obtained during the magnetic 
resonance image examination; at least one first processor for reconstructing the plurality of 

20 acquisition volumes to obtain a plurality of reconstructed volumes; and at least one second 
processor for registering the plurality of reconstructed volumes to obtain a plurality of 
registered volumes. The apparatus further comprises: means for processing the plurality of 
registered volumes to obtain statistical information; and means for outputting the statistical 
information. 

25 Moreover, the invention includes an article of manufacture as a computer-readable 

medium embodying code segments for a magnetic resonance image examination. The 
code segments comprise: code segments for controlling the magnetic resonance image 
examination to obtain a plurality of acquisition volumes; code segments for controlling at 
least one first processor to reconstruct the plurality of acquisition volumes to obtain a 

30 plurality of reconstructed volumes; and code segments for controlling at least one second 
processor to register the plurality of reconstructed volumes to obtain a plurality of 
registered volumes. 
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The apparatus of the invention includes a computer programmed with software to 
operate the computer in accordance with the invention. Examples of "computer" include: a 
general purpose computer; an interactive television; a hybrid combination of a general 
purpose computer and an interactive television; and any apparatus comprising a processor, 
5 memory, the capability to receive input, and the capability to generate output. 

The article of manufacture of the invention comprises a computer-readable medium 
embodying code segments to control a computer to perform the invention. Examples of a 
"computer-readable medium" include: a magnetic hard disk; a floppy disk; an optical disk, 
such as a CD-ROM or one using the DVD standard; a magnetic tape; a memory chip; a 
10 carrier wave used to carry computer-readable electronic data, such as those used in 

transmitting and receiving electronic mail or in accessing a network, such as the Internet or 
a local area network ("LAN"); and any storage device used for storing data accessible by a 
computer. Further, examples of "code segments" include: software; instructions; computer 
programs; or any means for controlling a computer. 



15 



BRIEF DESCRIPTION OF THE DRAWINGS 



The file of the patent contains at least one drawing executed in color. Copies of 
this patent with the color drawing(s) will be provided by the Patent and Trademark Office 
upon request and payment of the necessary fee. 
20 Embodiments of the invention are explained in greater detail by way of the 

drawings. 

Figure 1 illustrates a schematic of the hardware configuration of the invention. 

Figure 2 illustrates a flowchart of the software configuration of the invention. 

Figure 3 illustrates the processing configuration of real-time interactive fMRI using 
25 four processors. 

Figure 4. illustrates the processing configuration of real-time interactive fMRI 
using two processors. 

Figure 5 A shows plots for normalized standard error histograms of all vexed signal 
intensities for three BOLD fMRI experiments using the invention. 
30 Figure 5B shows plots for the degree of rotation detected and corrected using 

COCGV for the three BOLD fMRI experiments of Figure 5 A. 

-10- 
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Figure 6 shows images produced by the invention within 20 seconds after 
completion of a BOLD fMRI exam. 

Figure 7 shows images produced by the invention with 49 seconds after completion 
of a dynamic contrast fMRI exam. 

5 

DETAILED DESCRIPTION OF THE INVENTION 

The invention is described in terms of the hardware configuration, the software 
configuration, the processing configuration, an example, and the results from using the 
invention. 

10 Hardware Con figuration 

A preferred embodiment of the hardware configuration is illustrated schematically 
in Figure 1 . A clinical 1 .5 Tesla Echo Speed MR unit 1 (General Electric, Milwaukee, WI) 
is equipped with fast gradient hardware and a 2 gigabyte (GB) magnetic disk drive. A Bit3 
VME bus adapter card 2 (Bit3 Computer Corporation, Minneapolis, MN) was added to the 

15 MR unit 1 to enable access to the raw image data stored in the bulk access memory (BAM) 
3. An S-Bus Bit3 board 4 was also installed on a SPARC 20 Clone 5 (Axil Computer Inc., 
Concord, MA) and connected to the Bit3 card 2 on the MR unit 1. The SPARC 20 Clone 5 
is connected to a computational server 6 using a 10 Base T connection 7. The Bit3 
connections facilitated rapid transfer of raw image data sets from the MR unit 1 to the 

20 computational server 6 via the SPARC 20 Clone. The computational server 6 includes a 
four-processor DEC Alpha 4100 (Digital Equipment Corporation, Maynard, MA). In 
particular, the computational server is equipped with four 433 MHz computer processing 
units (CPUs), 1 GB of shared random access memory (RAM), and a 1 6 GB magnetic disk 
drive. 

25 The SPARC 20 Clone 5 plays two major roles. First, to synchronize the data input 

and output, the clone 5 serves as a data buffer between the MR unit 1 and the 
computational server 6. Second, the clone 5 performs byte-swapping of the data from the 
scanner in the MR unit 1 . The second role is necessary due to the incompatibility in the 
representation of the data between the scanner and the computational server. Further, with 

30 additional hardware (e.g., additional SPARC 20 Clones and Bit3 cards), multiple MR units 
can be connected to the computational server 6, which will allow for time-sharing of the 
hardware configuration. 
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The hardware configuration of the preferred embodiment is not a unique solution 
because it is designed to interface with a specific MR unit hardware and its output. As one 
of ordinary skill will recognize, analogous computer hardware configurations can be 
utilized with alternative clinical scanners. 

5 

Software Configuration 

In Figure 2, a flowchart illustrates the method carried out by the software of the 
invention. In block 8, raw image data is acquired by the MR unit 1 . The user interface 
preferably allows for the selection of the type of MR pulse sequence used (e.g., EPI or 
10 spiral imaging), the type of study being performed (e.g., BOLD dynamic contrast, spin 
tagging), number of volumes acquired, the matrix size, the paradigm definition, and the 
statistical thresholds. 

In block 9, the raw image data from the MR unit 1 is transferred to the RAM of the 
computational server 6 using appropriate software for transferring data. Moreover, the 
15 software preferably accepts raw image files directly from the MR unit 1 or from storage on 
the hard disk, tape, or other computer-readable medium. 

In block 1 0, the images are reconstructed by interpolating single-shot spiral data 
sets on to an orthogonal equidistant grid (27), followed by a two-dimensional Fast Fourier 
Transformation (2DFFT). The interpolation consists of a regridding algorithm with a 
20 Gaussian convolution window. 

In block 1 1 , the images are registered using a high speed optimized version of the 
correspondence of closest gradient voxels (COCGV) registration algorithm (28). This 
algorithm aligns areas of high intensity change based upon a rigid body transformation of 
the data sets. Of note, this particular regridding procedure allows for reconstruction of data 
25 collected with arbitrary k-space trajectories, including the trajectories characteristic of 
spiral scanning and EPI (24).. 

In block 12, statistical analysis of the images is performed. For BOLD fMRI 
studies, a normalization of the signal intensities across the volumes (VN) is performed. All 
volumes are normalized so that the mean intensity value of their relevant voxels is 
30 identical, which allows for comparison between studies. A baseline correction (BC) is 
applied on a voxel-by-voxel basis for all voxels within the subject using a third order 
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polynomial fit over the time series data (17). The "activation" statistical map is determined 
using an arbitrary user defined threshold for a student t-test. 

In block 13, the results of the fMRl exam are displayed for review by the 
investigator. All activation maps and graphs are displayed on a local video monitor in the 

5 MR operator control room ; or a hard copy is available using an interface to a color deskjet 
printer (Hewlett Packard, Palo Alto, CA). In addition, the results can be saved to a 
computer-readable medium or transferred over a network. 

Voxels determined to be statistically significant are superimposed in red onto the 
initial reconstructed fMRI spiral volume and displayed automatically. An s-map is also 

10 displayed and is defined as the standard error of the difference between the scan signal 
intensity on a voxel-by-voxel basis for all volumes in each epoch (i.e., on and rest 
conditions) in a box-car paradigm. Further, a histogram of the normalized standard error 
of all voxel signal intensity over time is plotted and displayed. The normalized standard 
error is the standard error of the difference between the two means (on vs. off) condition at 

15 each voxel divided by the mean signal intensity for these voxels (1 7, 1 8). 

Other types of information can also be determined and displayed using the 
invention. For example, for dynamic contrast studies, maps of rCBV, rCBF, and TTP are 
calculated by fitting a gamma variate function to the signal intensity time-course (1, 7, 8, 9, 
1 1, 29, 30). The fitting is performed for the voxels (50,000 to 150,000) within the subject. 

20 The TTP map is obtained by an algorithm searching for the occurrence of the maximum 
change in signal intensity over time, on a voxel-by-voxel basis. The arterial input function 
is obtained by searching all voxels within the subject for the greatest signal intensity 
change and then selecting the voxels with the earliest TTP signal change (8). The average 
gamma variate fit for the voxels meeting these two criteria was then used for the arterial 

25 input function. The deconvolution is based on singular value decomposition (29, 30). The 
rCBV, rCBF, and TTP maps are then displayed for review by the investigator. 

As another example, CBF images acquired by either arterial spin tagging (AST) or 
a multislice flow-sensitive alternating inversion recovery (FAIR) are processed in a similar 
manner to the BOLD fMRJ studies (1 2-16, 3 1). For these studies, a Tl map is crated 

30 during a separate acquisition, and these images are registered to the fMRI study. As with 

the BOLD fMRI studies, the percent change in CBF can be determined by taking the 

difference between the activation and rest states in response to a neurophysiological 
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stimulus. The CBF images and the Tl map are then displayed for review by the 
investigator. 

As will be apparent to one of ordinary skill, various other software routines can be 
implemented to accomplish the functions illustrated in Figure 2. 

5 

Processing Software Configuration 

Figure 3 illustrates a four parallel processor implementation of the invention using 
the hardware and software configurations described above. The software for the parallel 
processors was written in C++ with a TCL/TK interface or graphical interface. The 

10 scanner of the MR unit 1 acquires N volumes of raw image data, numbered V 1 , V2, V3, 
VN in Figure 3. The first acquired volume VI is transferred to the first processor PI of the 
computational server 6. The first processor PI reads and reconstructs the first volume VI, 
which is then used as the baseline volume for the registration of the remaining N-l 
acquisition volumes to correct for motion during the scan. 

15 The second acquired volume V2 is transferred to the second processor P2 of the 

computational server 6. The second processor P2 reads and reconstructs the second 
volume V2, and passes the reconstructed volume V2 to the third processor P3 of the . 
computational server 6. The third processor P3 registers the reconstructed volume V2 
against the reconstructed volume VI . 

20 The third acquired volume V3 is transferred to. the first processor PI of the 

computational server 6. The first processor PI reads and reconstructs the third volume V3, 
and passes the reconstructed volume V3 to the fourth processor P4 of the computational 
server 6. The fourth processor P4 registers the reconstructed volume V3 against the 
reconstructed volume VI. 

25 The fourth acquired volume V4 is transferred to the second processor P2 of the 

computational server 6. The second processor P2 reads and reconstructs the fourth volume 
V4, and passes the reconstructed volume V4 to the third processor P3 of the computational 
server 6. The third processor P3 registers the reconstructed volume V4 against the 
reconstructed volume VI . 

30 The above process continues for reading, reconstructing, and registering all N 

acquisition volumes. If N is an even integer, processor PI reads and reconstructs acquired 

volume VN-1, and processor 2 reads and reconstructs acquired volume V N . If N is an odd 
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integer, processor P2 reads and reconstructs acquired volume V N -i, and processor PI reads 
and reconstructs acquired volume V N . 

After all the N acquisition volumes have been read, reconstructed, and registered by 
the four processors P1-P4, the statistics of the fMRI examination are analyzed by the 
5 computational server. If one statistical analysis is to be performed, this statistical analysis 
can be performed by one or more of the processors. Moreover, if a plurality of statistical 
analyses are to be preformed, these statistical analyses can be performed by a combination 
of the processors. 

Finally, the resulting images, maps, and plots are displayed. The processors of the 

10 computational server are utilized as necessary to produce the displayable results. 

As one of ordinary skill will recognize, the paradigm illustrated in Figure 3 can be 
extended to any number of processors. For example, in Figure 4, the invention is 
illustrated using two processors. As another example, the paradigm can be used with 2*N 
processors, where N is the number of volumes acquired by the MR unit. For this example, 

15 each of the first N processors reads and reconstructs one of the N acquired volumes, and 
each of the second N processors registers one of the N reconstructed volumes. 

As will be apparent to one of ordinary skill the invention can be implemented 
using any number of processors and a wide variety of MR techniques. For example, the 
paradigm can be used with BOLD fMRI, dynamic contrast fMRI, the profusion imaging 

20 MRI technique, and cardiac MRI. Moreover, the invention can be used to combine 

diffusion imaging with fMRI data, or to combine spectroscopic imaging with fMRI data. 

Further, as will be apparent to one of ordinary skill, the invention is scalable as 
faster MR acquisition techniques become available and as larger data files are created. For 
example, in phased array acquisition, the raw image data generated by the multiple 

25 receivers of the phase array can be processed using multiple processors according to the 
invention. For instance, if four receivers are used, each of the data streams from each of 
the four receivers can be processed using four processors, as illustrated in Figure 3, and a 
total of 16 processors can be used to process all four of the data streams from the phased 
array. 

30 Moreover, the invention can be used as a stand alone post-processing unit. 

Additionally, the invention can be used with other compatible software packages for a 

more extensive statistical analysis of the fMRI data set. 
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Example 

Using the invention, a BOLD fMRI exam and dynamic contrast study were 
performed on the brain of a subject. In order to allow for the highest image acquisition rate 

5 within the MR scanner's gradient hardware constraints, a multislice single shot spiral scan 
acquisition technique was used for all fMRI modalities. Using the maximum gradient 
slew-rate of 120 T/m/s along with the maximum gradient amplitude of 22 mT m* 1 , the 
minimal acquisition window was 22 ms (18, 24) for a 64 x 64 matrix size and a 240 mm 
field of view (FOV), and the maximum MR scanner acquisition rate was 30 slices/second. 

10 Under the same conditions, the duration of the acquisition window for EPI scans was 35 
milliseconds. For both the BOLD fMRI studies and the dynamic contrast studies, 
multislice single shot gradient echo spiral images were obtained. The whole brain was 
covered in 2 seconds per volume with a nominal in-plane resolution = 3.75x3.75 mm 2 , 
TE=35 ms, TR=55.6 ms, number of slices = 36, slice thickness = 4 mm, and flip 

15 angle = 84_. The total acquisition time for the BOLD fMRI study was 4.13 minutes, and 
the total acquisition time for the dynamic contrast study was 2.67 minutes. Finger tapping 
activation paradigms (17) were used to evaluate the system performance. For the dynamic 
contrast studies, Gadopentetate Dimeglumine (Magnevist, Berlex Laboratories, Cedar 
Knolls, NJ) at 0. 1 to 0.2 millimole/kg was injected at 10 cc/sec using an MRI compatible 

20 mechanical injector (Spectris, Medrad, Inc., Pittsburgh, PA). To obtain baseline images, 
there was a scan delay of approximately 15-20 seconds prior to the injection. All clinical 
studies were performed under an approved Intramural Review Board Protocol at the 
National Institutes of Health. 

25 Results 

For the BOLD fMRI study, Table 1 summarizes the performance characteristics of 
the invention illustrated in Figure 3 of dividing the various tasks over four processors. 
Approximately 4 minutes and 8 seconds are required to acquire all 124 raw data volumes 
used for the BOLD fMRI study. During the acquisition time, the system processes the raw 
30 image files with essentially no delay. All 124 volumes are reconstructed and registered, 
and the statistical analysis is performed and displayed with a running time of 4 minutes and 
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28 seconds. Hence, with the invention, the information was displayed to the investigator 
20 seconds after the fMRI examination was completed. 



10 



TABLE 1 

Performance Characteristics of System 
Configuration for BOLD fMRI Examination 



36 slice volume acquisition time: 


2.0 sec/vol 


Raw data volume size: 


792000 bytes 


DEC read volumes (processors 1, 2): 


0.5 sec/vol 


Reconstruction (processors 1,2): 


1.5 sec/vol 


Registration (processors 3, 4): 


1.5 sec/vol 


VN, BC, and statistics (processors 1,4): 


10 sec. 


Display t-maps with images: 


7 sec. 


a) Total fMRI acquisition (124 volumes): 


4 min. 8 sec. 


b) Total processing time with BOLD fMRI 
data acquisition: 


4 min. 28 sec. 


Processing time (b-a) after MR data collection ends: 


20 seconds* 


c) Total processing time with BOLD fMRI 
data acquisition using a single DEC 
processor: 


15 min. 38 sec. 


Time ratio (c/b): 


3.5 



VME bus. 



The invention provides the ability to produce and display immediately the 
activation t-map, s-map, and normalized standard error histogram 20 seconds after the 
completion of the MR data acquisition. Further, the BOLD fMRI study results are reliable 
and reproducible, and are consistent with activation maps obtained with conventional 
15 processing strategies, which take hours of processing time (21 , 25). 

The total processing time reported in Table 1 must be compared to the several 
minutes of "down-time" usually experienced when performing fMRI studies. This delay is 
either due to waiting for the array processor to finish reconstructing the images or for 
writing the unprocessed raw data files to disk and subsequent transfer of the acquired data 
20 files (e.g., typically greater than 100 megabytes) to an off-line workstation. 

As an example of comparing the invention with the conventional approach. Table 1 
also includes the results from performing the same BOLD fMRI study, but using instead a 
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single DEC processor. Using the single processor, the computation time to process the 
large raw data BOLD fMRI files (i.e., 124 volumes of 64 x 64 by 36 slices) and produce t- 
maps was 15 minutes and 38 seconds. This included the time for disk retrieval of raw data 
using the single processor and the processing time using the software described above. In 

5 contrast, when all four processors on the computational server were used, the total 

processing time was 4 minutes and 29 seconds. Hence, using all four processors results in 
approximately a 3.5 fold increase in the speed of the data processing and analysis. 

Table 2 summarizes the performance characteristics of the invention illustrated in 
Figure 3 used with the dynamic contrast fMRI study to calculate rCBV, TTP, and rCBF 

10 maps. The initial steps in the process are exactly the same as for the BOLD fMRI study, 
except that fewer whole brain volumes (80) are acquired over the 2 minutes and 40 seconds 
acquisition time. After completion of the MR data acquisition, it takes an additional 49 
seconds using the invention to perform a voxel-by-voxel gamma variate fit of the time 
series data, search for the bolus TTP, determine the arterial input function used to calculate 

15 the rCBF maps, and display all three maps. Moreover, when the raw data dynamic contrast 
studies are processed from the disk using a single processor, the total processing time was 
10 minutes and 8 seconds. Hence, with a single processor, a delay of 8 minutes and 8 
seconds occurs, whereas with the invention, a delay of 49 seconds occurs. 

20 TABLE 2 

Performance Characteristics of System 
Configuration for Dynamic Contrast fMRI Examination 



36 slice volume acquisition time: 


2.0 sec/vol 


Raw data volume size: 


792000 bytes 


DEC read volumes (processors 1,2): 


0.5 sec/vol 


Reconstruction (processors 1,2): 


1.5 sec/vol 


Registration (processors 3, 4): 


1.5 sec/vol 


Calculation of rCBV, rCBF, and TTP maps (processors 1-4): 


29 sec. 


Display of rCBV, rCBF, and TTP maps: 


17 sec. 


a) Total fMRI acquisition (80 volumes): 


2 min. 40 sec. 


b) Total processing time with dynamic contrast fMRI data 
acquisition: 


3 min. 29 sec. 


Processing time (b-a) after MR data collection ends: 


49 seconds* 


c) Total processing time with dynamic contrast fMRI data 
acquisition using a single DEC processor: 
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lOmin. 8 sec. 


Time ratio (c/b): 


2.91 



VME bus. 

Using the invention, a BOLD fMRI exam was performed on a healthy control 
subject who was performing a sensorimotor activation by finger tapping at 2 Hz with the 

5 dominant hand. The data was processed using the invention illustrated in Figure 3. Figure 
6 shows images produced by the invention for the finger tapping. In Figure 6, activated 
voxels in red are superimposed on top of axial spiral images for the five motor regions (i.e., 
primary sensorimotor, lateral premotor region, parietal region, supplementary motor area, 
and ipsilateral cerebellum) associated with right-handed finger tapping. With the 

10 invention, the images were obtained within 20 seconds after acquiring the raw image data 
from the BOLD fMRI exam using a threshold of 4.8. 

With the invention, the normalized standard deviation histogram, along with a 
graph of patient motion which was obtained from the registration of the images, can be 
used as the criteria for rejecting a run or determining the necessity of repeating a given 

15 experimental run because of, for example, motion of the subject. In general, subject 
motion is difficult to address because, even when a registration finds the correct motion 
transformation, the ensuing interpolation reduces the high frequency components in the 
registered columns, which increases the standard deviation at each voxel position (32). 
To illustrate this aspect of the invention, the BOLD fMRI examination was 

20 repeated three times during the fMRI exam session. In the first session, the subject was 
asked to remain motionless. In the second session, the subject was asked to move during 
the exam to simulate an uncooperative patient. In particular, the subject was asked to 
begin the session with the subject's head tilted to one side and then to rotate it slowly to the 
other side during the course of the experiment. To demonstrate the utility of the interactive 

25 nature of the invention, the subject was again asked to remain motionless, and a third 
BOLD fMRI data set was subsequently collected and analyzed. 

Figure 5 A shows the normalized standard deviation histograms from the three 
sequential studies. As can be seen, the histogram from the second BOLD fMRI exam has a 
wider distribution and is interpreted as an inadequate exam precluding any further analysis 

30 of the data set. 
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Figure 5B is a plot of the rotational motion correction applied by the COCGV 
registration for the three BOLD fMRI experiments. All subject motion is represented as a 
single rotation magnitude and a single translation magnitude. The axis of rotation and 
translation are not displayed. 

Figure 7 shows rCB V, TTP, and rCBF maps produced by the invention for a 
dynamic contrast fMRI study. In Figure 7, each of the three rows shows six representative 
slices from a 36-slice whole brain dynamic contrast fMRI study performed on a healthy 
25-year old male control volunteer following administration of 0.1 mm/kg Gadopentetate 
Dimeglumine as a bolus at 10 cc/sec. using an MR-compatible injector. Row A shows the 
resulting rCBV maps, row B shows the resulting rCBF maps, and row C shows the 
resulting TTP maps. These maps were produced from 80 whole brain spiral images and, 
with the invention, were displayed within 49 seconds after completion of the dynamic 
contrast study. Further, approximately 99% of the voxels were successfully fitted by the 
gamma variate function. These maps are consistent with results obtained from a healthy 
control subject. 

The invention has been described in detail with respect to preferred embodiments, 
and it will now be apparent from the foregoing to those skilled in the art that changes and 
modifications may be made without departing from the invention in its broader aspects, 
and the invention, therefore, as defined in the appended claims is intended to cover all such 
changes and modifications as fall within the true spirit of the invention. 
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CLAIMS 



What is claimed is: 

1 . A method for outputting examination results from a magnetic resonance 
5 image examination, comprising the steps of: 

obtaining first and second acquisition volumes from a magnetic resonance image 
scanner during the magnetic resonance image examination; 

processing the acquisition volumes from the magnetic resonance image scanner to 
obtain the examination results, comprising the sub-steps of: 
10 reconstructing the first acquisition volume using a first processor to obtain a 

first reconstructed volume; 

reconstructing the second acquisition volume using a second processor to 
obtain a second reconstructed volume; and 

registering the second reconstructed volume using a third processor to 
15 obtain a second registered volume; and 

outputting the examination results. 

2. The method of claim 1 , further comprising the steps of: 

obtaining a third acquisition volume from the magnetic resonance image scanner 
20 during the magnetic resonance image examination; and 

wherein processing acquisition volumes further comprises the sub-steps of: 

reconstructing the third acquisition volume using the first processor to 
obtain a third reconstructed volume; and 

registering the third reconstructed volume using a fourth processor to obtain 
25 a third registered volume. 

3. The method of claim 2, further comprising the sub-steps of: 

obtaining a fourth acquisition volume from the magnetic resonance image scanner 
during the magnetic resonance image examination; and 
30 wherein processing acquisition volumes further comprises the sub-steps of: 

reconstructing the fourth acquisition volume using the second processor to 
obtain a fourth reconstructed volume; and 
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registering the fourth reconstructed volume using the third processor to 
obtain a fourth registered volume. 



5 



4. The method of claim 1 , further comprising the step of: 

processing the examination results to obtain statistical results of the magnetic 



resonance image examination using a combination of the first, second, and third 
processors. 

5. The method of claim 1 , wherein the magnetic resonance image examination 
10 comprises one of a blood oxygenation level dependent (BOLD) based functional magnetic 
resonance imaging (fMRI) examination, a dynamic contrast fMRI examination, a profusion 
imaging fMRI examination, a cardiac fMRI examination, diffusion imaging examination, 
and a spectroscopic imaging examination. 

15 6. The method of claim 1 , wherein outputting the examination results 

comprises at least one of displaying the examination results, saving the examination results 
to a computer-readable medium, transferring the examination results over a network, and 
printing the examination results. 

20 7. The method of claim 1 , wherein the examination results comprise at least 

one of a magnetic resonance image, an activation map, a relative cerebral blood volume 
(rCBV) map, a relative cerebral blood flow (rCBF) map, a time-to-peak (TTP) map, an 
activation s-map, an activation t-map, and an activation error histogram. 

25 8. A method for outputting examination results from a magnetic resonance 

image examination, comprising the steps of: 

obtaining first and second acquisition volumes from a magnetic resonance image 
scanner during the magnetic resonance image examination; 
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reconstructing the second acquisition volume using the first processor to 
obtain a second reconstructed volume; and 

registering the second reconstructed volume using a second processor to 
obtain a second registered volume; and 

outputting the examination results. 

9. An apparatus for outputting examination results from a magnetic resonance 
image examination, comprising: 

means for obtaining first and second acquisition volumes from a magnetic 
resonance image scanner during the magnetic resonance image examination; 

means for processing the acquisition volumes from the magnetic resonance image 
scanner to obtain the examination results, comprising: 

means for reconstructing the first acquisition volume using a first processor 
to obtain a first reconstructed volume; 

means for reconstructing the second acquisition volume using a second 
processor to obtain a second reconstructed volume; and 

means for registering the second reconstructed volume using a third 
processor to obtain a second registered volume; and 
means for outputting the examination results. 

1 0. An apparatus for processing data from a magnetic resonance image 
examination, comprising: 

a magnetic resonance image scanner for producing a plurality of acquisition 
volumes obtained during the magnetic resonance image examination; 

at least one first processor for reconstructing the plurality of acquisition volumes to 
obtain a plurality of reconstructed volumes; and 

at least one second processor for registering the plurality of reconstructed volumes 
to obtain a plurality of registered volumes. 

1 1 . The apparatus of claim 1 0, further comprising: 

means for processing the plurality of registered volumes to obtain statistical 
information; and 
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means for outputting the statistical information. 

12. A computer-readable medium embodying code segments for a magnetic 
resonance image examination, the code segments comprising: 

code segments for controlling the magnetic resonance image examination to obtain 
a plurality of acquisition volumes; 

code segments for controlling at least one first processor to reconstruct the plurality 
of acquisition volumes to obtain a plurality of reconstructed volumes; and 

code segments for controlling at least one second processor to register the plurality 
of reconstructed volumes to obtain a plurality of registered volumes. 
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